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Abstract 

Stochastic description of inflationary spacetimes emulates the growth of vac- 
uum fluctuations by an effective stochastic "noise field" which drives the dy- 
namics of the volume-smoothed inflaton. We investigate statistical properties 
of this field and find its correlator to be a function of distance measured in 
units of the smoothing length. Our results apply for a wide class of smoothing 
window functions and are different from previous calculations by Starobinsky 
and others who used a sharp momentum cutoff. We also discuss the appli- 
cability of some approximate noise descriptions to simulations of stochastic 
inflation. 
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I. INTRODUCTION 



Inflationary models explain structure formation as a consequence of vacuum fluctuations 
of a scalar field 0, the inflaton [|l| . Quantum fluctuations are also responsible for the structure 
of the universe at super- large scales and for the eternal character of inflation Because 
of exponentially rapid expansion of the spacetime, fluctuations of the inflaton field <fi on 
super- horizon scales effectively "freeze" a few Hubble times H~ l after they leave the horizon 
and can be treated as contributions to classical values of the inflaton field. After smoothing 
on super-horizon scales, the averaged field $ can be considered to be classical with an 
approximate equation of motion driven by the inflaton potential V ($) and by small-scale 
"noise" f , 



* = -g^(*)+e(x,t). a) 

Here the "noise field" £ (x, t) is a Gaussian random field characterized by its two-point 
function (£ (x, t) £ (x', t')). It is important to know the correct behavior of this correlator, 
in particular for numerical simulations of eternal inflation. 

One can approximate the noise field by using Eq. ([I]) with a free scalar field in de Sitter 
spacetime, even though the scalar field is only approximately massless and the spacetime is 
approximately de Sitter. The results will be applicable on scales smaller than the character- 
istic scale of the variation of the expansion rate H and that scale is typically exponentially 
large. In the pioneering work of Starobinsky 0] the correlator was evaluated in this manner 
as 



The smoothing scale for the inflaton field was R = [eHexp(-Ht)] 1 with e < 1. Eq. (|) 
shows a surprisingly slow decay of correlations at large distances. For comparison, the 
two-point function of the time derivatives of the un-smoothed field (^) (x, t) (j) (x', t')\ at 
large separations r = |x — x'| behaves as cx r~ 4 (here the angular brackets denote vacuum 
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expectation value rather than statistical average). One would not expect a smearing of the 
field operators </> (x, t) on scales R to have such an effect on correlations at distances r S> R. 

Our analysis shows that the origin of the unusual behavior of the correlator found by 
Starobinsky is the sharp momentum cutoff in his smoothing procedure. With a smooth 
cutoff, we recover the r -4 behavior independently of the cutoff window function and find 
that the time dependence of the noise correlator at large times is generically oc exp (—2Ht) 
instead of a sharp 5- function dependence of Eq. @. 

The correlator of noise is also important for simulations of the inflating spacetime using 
Eq. (HD- One of the simulation methods introduced in Ref. []|] consists of approximating 
Eq. (|]) by a finite difference equation with random sine waves playing the role of £, 

4> (t + At, x) - (t, x) = - — At + A sin (#e m nx + a) . (3) 

3x3 ^ ' 

Here n is a randomly directed unit vector, a is a random phase and A is an appropriately 
distributed random amplitude. We investigate the limits of validity of this approximation 
and show that for any finite At, the long-distance correlation of the effective noise is asymp- 
totically the same as that given by Eq. (|2[) which, as we argue, is unphysically large. Even 
with At —>■ 0, the sum of independent sine waves over one Hubble time H~ l produces an 
effective noise field with correlations oc r~ 2 . This method therefore can only be used in 
contexts that do not require precise simulation of the large-scale correlations of noise. In 
cases where such correlations are important, the noise field may be simulated as a Gaussian 
random field with known correlator; this method was employed in Ref. ||. 

The paper is organized as follows. The results regarding the noise correlator are presented 
in Section [TT], and the analysis of the sine wave approximation of Eq. (|3|) is in Section [[II]. The 
necessary details of the calculations are delegated to appendices. Appendix |A] describes the 
properties of window functions used for spatial smoothing, Appendix [FJ contains a deriva- 
tion of the noise correlator with arbitrary window, and Appendix |C| lists the correlators of 
unsmoothed field in de Sitter space for reference purposes. 
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II. SPATIAL AVERAGING AND NOISE 



Field fluctuations on super-horizon scales behave effectively as classical fluctuation modes 
with random amplitudes. This is conventionally described P,|T|-|IT| by averaging the field 
(f) in space over super-horizon scales and treating the resulting field as a classical (albeit 
stochastic) field $ satisfying an equation of motion with a stochastic term (or "noise field"). 
This term is a Gaussian random field which describes the effective character of the field 
fluctuations; it will be the focus of this section. 

We consider a free massless scalar field <fi in an inflating (de Sitter) spacetime with horizon 
size H^ 1 and the scale factor a (t) = exp (Ht). The equation of motion for the free field <fi is 

□0 = + 3#0--^A0 = O. (4) 
a [t) 



The field is quantized via the usual mode expansion (see Eqs. (|B1|) , ( p2|) of Appendix [BD, 

(x, t)= I (a^ k (t) e* kx + h.c) . (5) 

The averaging of the field <fi is performed by means of a suitable smoothing window W s (x; R) 
with a characteristic smoothing scale R, 

4> (x, t) = J rf 3 x> (x', t) W s (x - x'; R) . (6) 

Here, the physical smoothing scale is taken to be e _1 times larger than the horizon size, with 
e < 1. The corresponding comoving scale is 

The volume-averaged field has a mode expansion given by Eq. ( |A5| ) , 

r d 3 k 

4> (x, t) = / 3/2 w (kR) a k Vfc (t) e ikx + h.c, (8) 
J (2tt) 

where w (kR) is a suitable Fourier transform of the window profile W s (see Eqs. ([Al]), ( |A4| ) ) . 

The volume-averaged inflaton field is treated as a classical field $ satisfying Eq. (|IJ) 
which is a Langevin equation describing dynamics driven by an effective "noise field" £ (x, t) 



as well as by the effective potential V(<&). The noise field £ (x, t) in that equation can be 
heuristically defined as a stochastic field that corresponds to the quantum operator of the free 
field derivative </>, in the sense that any averages of £, such as the correlator (£ (x, t) £ (x', £')), 
are taken to be the same as the corresponding quantum expectation values of in the vacuum 
state (the standard Bunch-Davies vacuum). The effective noise field £ defined in this way 
is a Gaussian random field with zero mean, so the correlator (£ (x, i) £ (x',t')) completely 
describes its properties. The calculation of the noise correlator (£ (x, £)£ (x', £')) given in 
Appendix |B| is a straightforward computation of the corresponding expectation value of the 
quantum "noise operator" 0. 

The noise correlator generally depends on the particular window function W s (x; R) and 
on the parameter e. Since the smoothing window is merely a technical device in this ap- 
proach, we would expect to obtain results independent of the window W s (x; R) on scales 
^> R for a reasonably wide class of window shapes. It seems natural to require that the 
window function should be non-negative, spherically symmetric and depend on x only 
through the combination |x| /R. The last requirement makes W s a function of the form 
W s (x; R) = R~~ 3 W (|x| / R) with some suitable profile function W (q) which starts to de- 
crease rapidly at q ~ 1 (see Eqs. (|AT|)-([A2D). Additionally, we require that the window 
profile decays rapidly enough at large distances, as detailed below. 

It turns out, as shown in Appendix [B], that the asymptotics of the noise correlator at 
large temporal or spatial separations are independent of the choice of the window profile 
W (r) if it satisfies these conditions. The positivity assumption, W (r) > 0, may in fact be 
relaxed and substituted by the condition 




(9) 



For the integral in Eq. (^) to converge, we require that the window profile decays at large 
distances as W (r) ~ r -6 or faster. 



The noise correlator is a function of distance and of the smoothing scale e 



(£(x,t = 0)£(x' ) t)} = C(|x-x'|,t; e). 



(10) 
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While the precise shape of this function does depend on the window profile W(r), the 
large-distance asymptotic does not (cf. Eq. ( [B10| )): 



C(|x-x'|,t; e) oce- 2m (eH|x-x / |)~* + 0(|x-x'|~ °j . (11) 

This behavior agrees with the r -4 asymptotic behavior of the field derivatives in de Sitter 
space (see Eq. flC6f) ). However, Eq. ( |TT1) is in disagreement with the result of Starobinsky || 



where the correlator was found to have a much higher correlation of noise at large distances 
(Eq. ©)• 

The discrepancy between Eqs. (0) and ( |TTD is due to different choices of the smoothing 
window profiles. Ref. H used a sharp step-function cutoff in Fourier space, 

w (kR) = 9 (1 - kR) (12) 

which corresponds to the real-space smoothing window 

sin x — xcosx , . 

w ® = w ■ (13) 



The smoothing window of Eq. ( fIB| ) is not everywhere positive and decays too slowly at 
large distances to satisfy Eq. (|9[). The latter condition is equivalent to the requirement that 
w (kR) be a sufficiently smooth function (at least twice differentiable). If we allow windows 
W (x) with discontinuous Fourier transforms w (kR), then different choices of window would 
lead to a wide range of asymptotic behaviors of the correlator. On the other hand, as long 
as Eq. (P) is satisfied, the asymptotic Eq. ([TT|) is independent of the choice of the window 
shape w (kR) or equivalently W (x). Therefore we conclude that the long-distance behavior 
of the correlator obtained in Ref. |4]] is an artefact of the sharp mode cutoff. 

In the limit of small e we find that the noise correlator is a function only of p = eH |x — x'| 
(the distance measured in smoothing scale units) and t (see Eq. (|B12j )). It means that the 
correlator has a scaling property, 

C (|x - x'| , t; eA) = C (A |x - x'| ,t;e) + (e 2 ) . (14) 
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For illustration we give some specific results with a Gaussian smoothing window. The 
correlator profile in that case is given by Eq. ( |B13| ). The shape of the correlator is illustrated 
in Fig. 0. 



III. ON SIMULATIONS OF THE STOCHASTIC DYNAMICS 

Simulations of inflating spacetimes use a discretized version of Eq. (|4]). One then needs 
to represent the effective noise field £ by some random process. An approximation used in 
Refs. employed random sine waves in lieu of £, 

|(x,t) = Asm(Ha(t)nx + a) . (15) 

Here n is a randomly directed unit vector, a is a random phase and A is an appropriately 
distributed random amplitude. The equation of motion was discretized using the slow roll 
approximation Eq. ([!]), yielding 

4>{t + At, x) - (t, x) = - V '^ At + I (x, t) . (16) 

The time interval At was chosen to be much smaller than the Hubble time if -1 . 

The random field £ (x, t) defined by Eq. flT5] ) is not actually Gaussian, although the one- 
point distributions of £ (x, t) at any fixed point x are Gaussian if the random amplitude A 
is drawn from the standard \ 2 distribution. We can compute the correlator by averaging 
over A, a and n and find, up to normalization, 



(e(x,t)e(x',t'))oc 



sin H 


a (t) x — a (t') x'| 


H\ 


a 


(t)x- a(t')x'| 





(17) 



Note that the correlator at t = t' is of the same form as Eq. (|^), although there is no 
sharp time dependence. The field £ (x, t) exhibits the same large correlations ~ r _1 at large 
distances as the noise field computed with the sharp mode cutoff. 

One could hope to improve the approximation by adding many small timesteps At — » 
in Eq. (|HJ) while keeping the total time increment constant, NAt ~ if -1 . The net effect is 
that of adding many independent sine waves together, 
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1 N ~ l _, I rt+NAt _ 

£ (x, f ; iVAt) = - £ £ (x, t + nAt) « — jf £ (x, r) dr. (18) 

The correlator of the resulting random field £ (x, £; NAt) is 

£(x,t; ATAt)|(x',t'; NAt)) oc — (19) 

' H z \a (i) x — a (v) x'| 

(The numerator in the last equation is an oscillating function of order 1, the explicit form of 
which will not be important.) The field £ (x, t; NAt) is an approximation to an increment 
of the smoothed field A(p = (f) (x, t + NAt) — (x, t) , and the correlator in Eq. ( [T9|) with 
iVAt ~ if -1 is to be compared with the corresponding correlator of increments, which under 
the same assumptions as Eq. ( ]TT|) can be shown to behave at long distances as 

X A0 (x, t) A0 (x\ t')) oc |x - x'f 4 . (20) 



We find that the unphysically large correlations at large distances do not disappear even in 
the limit of At -> 0. 

We conclude therefore that the approximation method of Eq. (|l"5|) can only be used in 
contexts that do not require precise simulation of the large-scale correlations of noise. 
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APPENDIX A: PROPERTIES OF WINDOW FUNCTION 



We first investigate properties of the class of window profiles we will be using for spatial 
averages of (p. Spherical symmetry and scaling with R forces the smoothing window profile 
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W s (x; R) to depend on x only through the combination |x| /R, and from the normalization 
/ W s (x; R) d 3 x = 1 it follows that 



W a (x.;R) = Rr 6 W\^j (Al) 
with a suitable dimensionless function W. The normalization condition becomes 

An Jw (q) q 2 dq = 1. (A2) 

By construction, the window profile W (q) should fall off at q ~ 1 which corresponds to 
distances x of order H~ l . From Eq. (|A2|) it follows that W (q) must decay at least as q~ A 
at large q; however, we shall see below that the oc q~ 4 or even oc q~~ 5 decay of the window 
function introduces too much correlation between far-away points, and so we shall assume 
that W (q) falls off at least as g~ 6 or faster at large q. 

Smoothing is more conveniently performed directly in the Fourier domain, where it cor- 
responds to suppressing high-frequency modes of the field. The modes are attenuated by 
the Fourier image of the window function. From the form (WT) of the window function it 



follows that its Fourier image is a real function of kR (t) = p: 

e- ikx ^- 3 W (|) d 3 x = w (kR) , (A3) 
where 

w (p) = Ati r S -^-W (q) qdq. (A4) 
Jo p 

Since the profile W (q) starts to decay at q ~ 1, its Fourier image w (p) becomes negligible 
at p ^> 1. The mode expansion of the smoothed field is 

<p (x, t) = / w (kR) tfj k (t) a k e ikx + h.c. (A5) 

For instance, a Gaussian smoothing window in real space is 



W s (x; R) = exp 

(27r) 3/2 i? 3 



|x i 
2W- 



(A6) 



and in the Fourier domain 



k 2 R 2 ' 



w (kR) = exp 



exp 



k i] 



(A7) 



2 J \ 2e 2 J ' 

Restrictions on the window profile W (q) lead to conditions on its Fourier image w (p) . 
Normalization of the window profile means that w (0) = 1, while the fast decay of W (q) at 
large q makes w (p) a differentiable function, 



d n w (p) 



dp n 



An 



d n sin pq 



dp n p 



qW (q) dq 



<oo ifp ^ 0, \W(q)\ < Cq- n ~ 2 at q -> oo, 



(A8) 



d n w (p) 



dp r > 



<oo if \W{q)\ < Cg" atg -> oo, 



(A9) 



p=0 



where C is a suitable constant. The odd derivatives of w{p) at p = generally vanish: 



r/ 2 " • ' //• (0) 
— = An 

dp 2n+l 



f 

JO 



d 2n+l 



sin pq 



qW (q) dq = 0. 



(A10) 



p=0 



dp 2n+1 p 

(The latter condition is satisfied only if W (q) falls off faster than q- 2n ~ A at large q because 

only in that case it is legitimate to take the limit of p = within the integral.) We find that 

if W (q) is decaying faster than q~ A viz. q~ 5 then w' (0) = and w" (0) is finite. We shall 

see shortly that the first two derivatives of w(p) at p = and smoothness of w (p) at p > 

determine the asymptotic behavior of the correlators for the smoothed field. 

A consequence of the positivity of W (q) is the non-vanishing of the even derivatives of 

w (p) at p = (provided they exist), for instance 

47r r 00 



w "(0) = - q 4 W{q)dq<0. 
3 Jo 



(All) 



If W (q) falls off exponentially at q — > oo, the function w (p) can be expanded at small p 



as 



, w"(0) 2 ti/ 4 )(0) 4 
1 + — + — -j-^p 4 + 



Air r°° 

w (2n) (Q) = / g 2n + 2^ ( } (M2) 

2n + 1 Jo 



* ' 2! ^ ' 4! - ' -' " - ^ ' 2n + l7o 
The derivatives (0) depend on the window function but are generically of order 1 since 

W (q) falls off at q ~ 1 by construction. Our results depend only on the first two terms of 

this series and on the differentiability of w (p) at all p, so the condition W (q) > can be 

dropped as long as w" (0) < 0. 
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APPENDIX B: CORRELATOR OF NOISE 



Here we derive the correlators of the effective noise field £ (x, t) for an arbitrary smoothing 
window. 

The field <fi is quantized using the mode expansion 

d 3 k 



<t> (x, t) 



(27T) 



3/2 



(flk^fc (t) e ikx + h.c.) . 



(Bl) 



The mode functions are 



iZe-tf"? / 1 



(B2) 

where the conformal time rj is defined as 77 = — H~ l e~ Ht , with <9t = —Hrjd^. The mode 
functions ^ (r/) satisfy the equation of motion, i.e. 

d 2 .0 



fa (v) = 0. 



(B3) 



We have defined the noise field £ (x, t) through the time derivative of the averaged field 
which has a mode expansion 



0(x,t) 



d 3 k 



(27T) 



3/2 



v k (v) ake + /i.e. 



(B4) 



where 



v k ( V ) = Hj t [w (kR) ip k (77)] = H \-kRw' (kR) ^ fc (77) + w (kR) 4 (v) 



(B5) 



In the limit of e < 1 we may disregard the second term in the square brackets. The noise 
correlator is then 

H%r] 2 f 00 



(f (xi,7?i)£(x 2 ,772)) 



, , , dk sin fcr h (k) 
An 2 re 2 Jo 



h (k) = (1 + zA;77i) (1 — z&77 2 ) e 



iKm-vi) w > [ w > 



kr]2 



(B6) 



(B7) 



From Eq. ( |B(j| ) we can obtain the asymptotic form of the correlator at large effective 
distances p ^> 1. First we note that Eq. ( p6|) is a sine transform of a certain (dimensionless) 
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function which we denoted h (k). It is known that the asymptotic of a Fourier transform is 
determined by the degree of smoothness of the transformed function. In the particular case 
of the sine transform, repeated partial integration gives 
r°° h (0) h" (Q) 1 r°° 

/ h{k)smrkdk = ^--4^ + ^ dk cos kr ti" (k) . (B8) 
jo r r 6 H Jo 

If h'" (0) is finite, the last integral is O (r -5 ) and Eq. (|B8|) may be used as the asymptotic 

at large r. Therefore we need to compute the second derivative at k = of the relevant 

function h (k) in our case, and we can use the expansion (|A12|) for the window function 



w (p). We also find that the asymptotic expansion of Eq. ( p8| ) is valid, i.e. h'" (0) is finite, if 

«j'(0) = 0, |tw"(0)|<oo. (B9) 

As we have seen above, these conditions are satisfied if for example |VT(g)| < Cq~ & at 
q — > oo. The asymptotic form of Eq. (|B6|) at large r is then 



(X!, I,!) £ (x 2 , V2 )) = - ^JJlf W (0) | 2 + O (r- 6 ) . (BIO) 



This expression is very similar to Eq. (|C6|) for the unsmoothed correlator of field derivatives. 
We note that the asymptotic Eq. ( B10|) is essentially independent of the shape of the window 
function, since the value \w" (0)| as indicated by Eq. ( |A11| ) has the meaning of the window- 
averaged squared distance and must be of order 1 because the window profile W (q) starts 
to decay at q ~ 1 by construction. 

We can obtain a simpler expression for the correlator in the limit when the smoothing 
parameter e is small while the product eHr is finite. A rescaling r — > eHr = p and the 
corresponding change of variable k = eHu simplify Eq. ( |B6[ ) because we can omit terms of 
order e and smaller; in particular, the product of mode functions is simplified to 

V* (Vi) 1>k fa) = (l + O (e 2 )) • (Bll) 

The leading term in the correlator, expressed through k and p, becomes 

l£ (x x , 7&) £ (x 2 , V2 )) = | o °° dn sin np w' {-H Vi k) w' {-Ht] 2 k) + O (e 2 ) . (B12) 
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We find that in the limit of small e but finite eHr the correlator as function of the "effective 
distance" p and the time difference (expressed by 772/771) becomes independent of e. 

Eq. ( |B12| ) allows us to compute the correlator at all distances in the limit of small e. For 
a Gaussian smoothing window, w (p) = exp (—p 2 /2), we obtain in the limit of small e 



(f (xi,r/iK(x 2 , 77 2 )) 



exp 



k 2 sin Kpdn 



(if%r/ 2 ) 2 z/ 4 



v 1 



iv \ 

VV3j exp 



V" 



(B13) 



4vr 2 p 4 

where we introduced a dimensionless quantity 

v = P (B14) 

^H 2 ( V 2 + r/1) 

The shape of this function for 771 = r/ 2 is shown in Fig. 1. The leading term of the expression in 
brackets in Eq. ( |B13[ ) at large v is (— 2z/~ 4 ), and since for the Gaussian window w" (0) = — 1, 
we recover Eq. (B10| ). The value of the correlator at coincident points (p = 0) as function 
of time separation is (cf. Eq. (|G4|)) 



(£ (0,770 £(0,772)) 



H 4 



(B15) 



2vr 2 (^2 + ^2)2 g7r 2 cogh 2 HAf ■ 

One can also obtain the leading asymptotics of the unequal-time correlator at large time 
separations. We again start with Eq. (|B~6|) and assume that the time separation is much 
greater than the Hubble time, f]2/i]i = a 1 ^ 1- F° r simplicity we can take Hrji = —1. We 
use Eq. ( |A12| ) for w (a _1 /c) at small a _1 K (since the integration is effectively performed over 
a fixed finite range of k) to obtain 



(Z (xi,77i)C(x 2 ,77 2 )) 



H 2 w" (0) 



47r 2 e 3 a 2 Hr 



k 



(Ik k sin krw' I ^ ) e ik/H [ I + 1— ) + 0{ <r 



k_" 
H 



(B16) 



The integral in Eq. (|B16| ) is clearly time-independent. Therefore the correlator decays as 
or 2 = exp (—2Ht) with time separation at any fixed distance. 



13 



APPENDIX C: CORRELATORS OF UNSMOOTHED FIELDS 



The correlator of time derivatives of the field at equal times is 



(x,t)0(x',t) 



1 f , , 2 ., , ,s „ 1 r°° sin At • . . 2 



(27T) 



e ik(x-x') rf 3 k 



27T 2 jo r 



i\) k (t) kdk, (CI) 



where r = |x — x'| is the distance between points. (In this Appendix, angular brackets 
denote quantum expectation values.) The derivatives of the mode function are 



dipk (v) ikHrj 



which gives (after a regularization of the integral) at t = (ifr/ = — 1 

1 a 00 sin At 



(x,t)0(x',t) 



47T 2 JO 



\ (P roo \ 

-k 2 dk = —^-r^ / sinkrdk = - . 
47r 2 dr 2 7o 2w 2 r 4 



The unequal-time correlator at coincident points (r = 0) is found as 

x2 



0(0,^)0(0,^) 



1 



2tt 2 Jo 



(F 2 r/ir/ 2 )- 



4tt 2 Jo 



(C2) 



(C3) 



k 3 e -ik( V2 - m ) dk 



An 2 d V iJo 4tt 2 ( Vl - r/ 2 ) 4 32tt 2 \ 2 



sinn - 



Finally, we consider the general correlator at arbitrary points and times: 

1 



(xi,ti)0(x 2 ,t 2 ) 



2tt 2 



Q1T1 fCf 

^ (t)r k (o) — kdk 



(C4) 



H 4 



(V1V2) 



2 3 (7/1 - r] 2 f + r 2 



3 • 



2n ' ' ((V1-V2)'- 
As expected, it diverges on the lightcone r = \rji — 77 2 |. The asymptotic form of Eq. 



(C5) 



51) at 



large distances r is 



(xi,ti) 0(x 2 ,t 2 



# 4 fa^)' 
2ttV 4 



+ O (r~ 6 ) 



(C6) 
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FIGURES 
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FIG. 1. Correlator C (p,t) of the effective noise field as function of p computed at equal times 
(t = 0). The distance was measured in smoothing scale units (p = eHr). 
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